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Abstract 

We present here exact results for a one-dimensional gas, or fluid, of 
hard-sphere particles with attractive boundaries. The particles, which 
can exchange with a bulk reservoir, mediate an interaction between 
the boundaries. A two-dimensional lattice of such one-dimensional gas 
'columns' represents a discrete approximation of a three-dimensional 
gas of particles between two surfaces. The effective particle-mediated 
interaction potential of the boundaries, or surfaces, is calculated from 
the grand-canonical partition function of the one-dimensional gas of 
particles, which is an extension of the well-studied Tonks gas The 
effective interaction potential exhibits two minima. The first minimum 
at boundary contact reflects depletion interactions, while the second 
minimum at separations close to the particle diameter results from 
a single adsorbed particle that crosslinks the two boundaries. The 
second minimum is the global minimum for sufficiently large binding 
energies of the particles. Interestingly, the effective adhesion energy 
corresponding to this minimum is maximal at intermediate concentra- 
tions of the particles. 

1 Introduction 

The interactions of surfaces are often affected by nanoparticles or macro- 
molecules in the surrounding medium. Non-adhesive particles cause attrac- 
tive depletion interactions between the surfaces, since the excluded volume 
of the molecules depends on the surface separation [21 El S] - Adhesive parti- 
cles, on the other hand, can directly bind two surfaces together if the surface 
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separation is close to the particle diameter [SI El El- In a recent letter [S], 
we have presented a general, statistical-mechanical model for two surfaces in 
contact with adhesive particles. In this model, the space between the sur- 
faces is discretized into columns of the same diameter d as the particles. The 
approximation implied by this discretization is valid for small bulk volume 
fractions of the particles, since three-dimensional packing effects relevant at 
larger volume fractions are neglected. For short-ranged particle-surface in- 
teractions, the gas of particles between the surfaces is as dilute as in the 
bulk for large surface separations, except for the single adsorption layers of 
particles at the surfaces. 

In this article, we present an exact solution of the one-dimensional gas of 
hard-sphere particles in a single column between two 'surfaces'. Our aim here 
is two-fold. First, the exact solution presented here corroborates our previ- 
ous, approximate solution for this one- dimensional gas obtained from a virial 
expansion in the particle concentration [8]. Second, the exactly solvable, 
one-dimensional model considered here is a simple toy model to study the 
interplay of surface adhesion and particle adsorption. Exactly solvable, one- 
dimensional models have played an important role in statistical mechanics 
[HI [10]. One example is the Kac-Baker model [UJ [121 H2] , which has shed light 
on the statistical origin of phase transitions of the classical van der Waals 
type. More recent examples are models for one- dimensional interfaces, or 
strings, which have revealed the relevance of entropy and steric interactions 
in membrane unbinding and wetting transitions [Til E2> [IS] • Other examples 
are the Tonks model [1] and its various generalizations [T71 [181 HH1 EO] , which 
have influenced our understanding of the relations between short-ranged par- 
ticle interactions, thermodynamics, and statistical correlations in simple flu- 
ids. The Tonks model has been exploited also in soft-matter physics to 
investigate structures of confined fluids [2D [22], depletion phenomena in two- 
component mixtures [23] , thermal properties of columnar liquid crystals [23] 
and the phase behavior of polydisperse wormlike micelles [25] . A recent bio- 
physical modification of the Tonks model addresses the wrapping of DNA 
around histone proteins [2S] . The model considered here is a novel extension 
of the Tonks model. 

In our model, a one-dimensional gas of hard-sphere particles is attracted 
to the system boundaries, or 'surfaces', by short-ranged interactions. We 
calculate the effective, particle-mediated interaction potential between the 
surfaces, V, by explicit integration over the particles' degrees of freedom in 
the partition function. The potential V is a function of the surface separa- 
tion I and exhibits a minimum at surface contact, which reflects depletion 
interactions, and a second minimum at separations close to the diameter of 
the adhesive particles. The effective, particle-mediated adhesion energy of 
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Figure 1: a) A one-dimensional gas, or fluid, of hard-sphere particles in a 
column of the same diameter d as the particles. The interaction between 
the particles and the boundaries, or 'surfaces', is described by a square- well 
potential of depth U and range l r < d/2. A particle thus gains the binding 
energy U if its center is located at a distance smaller than l r + d/2 from one of 
the surfaces. The length I of the column corresponds to the separation of the 
surfaces. We consider the grand-canonical ensemble in which the particles 
in the column exchange with a bulk reservoir. - b) For small bulk volume 
fractions, a two-dimensional lattice of such columns represents a discrete 
approximation of a three-dimensional gas of particles between two surfaces 
[8]. Since the particle-surface interactions are short-ranged in our model, 
the particle gas between the surfaces is as dilute as in the bulk for large 
surfaces separations, except for the adsorption layer of particles at each of 
the surfaces. 



the surfaces, W, can be determined from the interaction potential V. The 
adhesion energy is the minimal work that has to be performed to bring the 
surfaces apart from the equilibrium state corresponding to the deepest well of 
the potential V(£). Interestingly, the adhesion energy W attains a maximum 
value at an optimal particle concentration in the bulk, and is considerably 
smaller both for lower and higher particle bulk concentrations. 

This article is organized as follows. In section 2, we introduce our model 
and define the thermodynamic quantities of interest. In section 3, we cal- 
culate the particle-mediated interaction potential V(£) of the surfaces. The 
global minimum of this potential is determined in section 4, and the effective 
adhesion energy of the surfaces in section 5. In section 6, we show that the 
interaction potential V(£) exhibits a barrier at surface separations slightly 
larger than the particle diameter, because a particle bound to one of the sur- 
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faces 'blocks' the binding of a second particle to the apposing surface. The 
particle binding probability is calculated and analyzed in section 7. 



2 Model and definitions 

We consider a one-dimensional gas of particles with attractive boundaries, 
see figure [TJ The particles are modeled as hard spheres, and the attractive in- 
teraction between the particles and the boundaries, or 'surfaces', is described 
by a square- well potential with depth U and range l r . The length I of the 
gas 'column' corresponds to the separation of the surfaces and the width of 
the column is chosen to be equal to the particle diameter d. The particles in 
the column exchange with a bulk reservoir of particles. 

The position of the center of mass of particle k is denoted by x k , and its 
momentum by p k . For the system of n hard particles confined in the column 
of length I > nd, one has d/2 < x\, xi<X2 — d, x 2 <X3~d, ... x n < 
I — d/2. We assume that the 1-st and n-th particle interact with the surfaces, 
i.e. with the bases of the columns, via the square-well potential 

v n {x k } = -ue \^ + ir- x?j - ue(x n -e + i r + ^), (1) 

where U > and l r > are the potential depth and range, respectively. We 
also assume that l r < d/2. Here and below, denotes the Heaviside step 
function with Q(x) = 1 for x > and = for x < 0. The configuration 
energy for the system of n particles in the column is 

n 2 

n n {x k , Pk } = v n {x k } + J2Y- (2) 

k=l 

and the corresponding canonical partition function can be written as 

i rl— f rx n -d rx?,~d rX2—d 

Zn = — dx n dx n ^ ... dx 2 d Xl e- V ^T (3) 

An J(n-l)d J(n-I)d J Id Jld 

after integration over the momenta of the particles, see, e.g., [H[2TJ[26]. Here, 
A = h/ [2-KmT) 1 / 2 is the thermal de Broglie wavelength, and T denotes the 
temperature times the Boltzmann constant. In other words, T is the basic 
energy scale. 

Since the particles can exchange with the bulk solution, the number n 
of particles in the column is not constant. Such a system is described by 
the grand-canonical ensemble in which the temperature, the column length 
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£, and the particle chemical potential \i are fixed. The corresponding grand- 
canonical partition function 



ie/d\ 

Z = l+J2 Z ^' T ( 4 ) 

n=l 

is a sum of a finite number of elements, where [£/d\ denotes the largest 
integer less than or equal Ijd. The upper limit of the sum on the right hand 
side of equation (j3j) is the largest number of hard particles with diameter d 
that can be accommodated in a column of length £. The partition function 
Z given by equation (@| determines the grand potential 

F gc = -TlnZ, (5) 

the bulk density of the grand potential 



/ gc = Jim *±* (6) 



and, hence, the surface contribution to the grand potential Fgc = F gc — 
f gc £/d. The effective interaction potential of the surfaces 0, 

p( s ) f? (j f p 

y _ r S c _ r gc a Jgc<- /_x 

d? d 3 { ' 

is defined as the density of the surface contribution to the grand potential 
F gc . For consistency with our previous model for particle- mediated surface 
interactions [8], the column bases are chosen here to be squares of side length 
d. Thus the d 2 in the denominator of equation (j7j) is the column base area. 
The surface potential V defined by equation (J^J) is the main quantity of 
interest here and will be determined in the next section. 



3 Effective surface interaction potential 

Equations (j4j) - ([7|) imply that 

exp (-V^ ) = ( 1 + g exp (Ipj 



1 The surface interaction potential V was called the 'effective adhesion potential' in 
reference [5]. 
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To determine the surface interaction potential V, we thus have to calculate 
the canonical partition function Z n and the bulk density of the grand po- 
tential, / gc , defined in equation Q). The n-particle partition function Z n 
is defined in equation ([3]). The change of variables yu = Xk — (k — |) d in 
equation ([3]), with k = 1,2 ... ,n, leads to 

i fl-nd />y n t-yz ryi 

Z n = i- n dyne ue(y n -t+n d+lr)/ T d^.../ dy J di/xe^-*)/* 

A Jo Jo 7o Jo 

(9) 

where [7 is the binding energy of the particles. The integral (Q can be 
evaluated, see [S], and after some computation we arrive at 

Z n = (e u ' T - l) 2 $ n (£ - 2l r ) - 2e u / T (e u ' T - l) $ n (i - l r ) + e 2U ' T $ n (I) , 

(10) 

with 

Mk) = i(^^Ve(Jo-nd) (11) 



n! V A 

for any length Z . For [/ = 0, equation ffTU]) . reduces to the partition function 

2„ = -^(^-nd) n e(^-nd) (12) 

of the classical Tonks gas P . The grand potential density f gc can be derived 
from this exact result as shown in the following subsection. 

3.1 Thermodynamic potentials in the bulk 

The canonical and grand-canonical ensembles are equivalent in the thermo- 
dynamic limit, i.e. for infinite surface separation I. In this limit, the grand 
potential density / gc therefore can be obtained from the canonical poten- 
tial density via Legendre transformation. First, we define the canonical free 
energy F ca = — T In Z n and the free energy density in the bulk, 

/ ca = lim-^— , (13) 



where lim^ denotes the thermodynamic limit in which both the column 
length i and the particle number n go to infinity while the particle 'volume 
fraction' 

nd , . 

= — (14) 

remains constant. The particle volume fraction in one dimension defined by 
equation f[T4|) attains values < <fi < 1, and could also be called a 'length 
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fraction'. In the one-dimensional model considered here, close packing cor- 
responds to (f) = 1. 

The free energy density / ca defined by equation (fT3|) is an intensive quan- 
tity in the thermodynamic limit and, therefore, does not depend on the 
boundary conditions. In particular, the free energy density / ca and its deriva- 
tives do not depend on the boundary potential (PJ, which is characterized by 
the binding energy U and range l r . With the exact expression for the canon- 
ical partition function Z n given in equation ffTOj) and the Sterling formula 



we get 



ln(n!) ~ n In n — n + — In (27rn) , (15) 



/ca = T(f) 



•1 + ln 



A 



(16) 



d 1 - 

As expected, the particle-surface interactions characterized by the binding 
energy U and range l r do not affect the free energy density / ca in the ther- 
modynamic limit. From the canonical free energy density / ca , we obtain the 
chemical potential 

" = {W) T (17) 

for the particles in the bulk. Equations (1161) and (IT7|) lead to 
which can be rewritten as 

Finally, the grand potential density / gc follows from the Legendre transfor- 
mation 

/gc = /ca-/i0. (20) 

Equations (fTBl). (JIg|) and (j20l lead to 



= (21) 

with the dependence between the chemical potential \x and the particle bulk 
volume fraction given in equation (ITS]) . Basic thermodynamics implies 
that the gas pressure in the bulk is p = — f gc /d 3 . This leads to the pressure 
p = T<f)/[(1 — <f))d 3 ] here, which is the correct equation of state for the Tonks 
gas pp. 
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3.2 Potential profile 



By combining equations (1191) and (f2~Tj) . the effective, particle- mediated 
interaction potential V of the surfaces can be expressed as a function of the 
particle bulk volume fraction and the separation i between the surfaces: 



T 



\tld\ 

n=l 



A 



d 1 - 



cxp 



n <p 



+ 



T 



d 3 1 - 



(22) 



With the exact expression for the canonical partition function Z n given in 
equations (TTOT) and (fTTi) . the potential V can be evaluated numerically for any 
finite surface separation £, see figures [2] and [3j For large bulk volume fractions 
<f> (dashed curve in figure[2]), the potential V exhibits oscillations up to surface 
separations I of the order of several particle diameters before approaching a 
constant, asymptotic value Voq for large surface separations. The oscillations 
are related to successive layers of particles formed in the space between the 
surfaces. For small bulk volume fractions 0, in contrast, the interaction 
potential V attains an approximately constant value for surface separations 
£ > 2(d + l r ), see solid curve in figure [2J 



3.3 Potential asymptote 

Since the interactions between the particles and the surfaces are short-ranged, 
the potential V has a horizontal asymptote, i.e. V(£) approaches a constant 
value Voo for large surface separation ft, see figures H] and El In this subsection, 
we calculate the position of the asymptote. To simplify the notation, we first 
introduce the auxiliary variable 



C = (23) 



and the function 



ft h A L ^.V ( t-nd a fi-i Y'fe-k 



n 



(24) 

defined for an arbitrary length / - By combining equations ([HD, (fIU|) . (ITS]) 
and (12 ip . we then express the potential V as 

e -v(e )d yT = e - W d + {e u/T _ 1)2 g ^ ^ c 
_ 2e ^(e^ - 1) g (i, l j, Cj + e 2VIT g 0, . (25) 
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Figure 2: Rescaled surface interaction potential Vd 2 /T, given by equation 
f[22]) . as a function of the rescaled surface separation tjd where d is the par- 
ticle diameter and T denotes the temperature in energy units. The particle 
binding energy here is U = 5T, the binding range is l r = 0.3 d, and the 
particle bulk volume fraction is = 0.05 (solid line) and = 0.5 (dashed 
line), respectively. For large surface separations £ d, the potential V(£) 
attains an approximately constant value V^. According to equation (1281) . 
the asymptotic values are ~ —6. 6446 Tjd 2 for = 0.5 (dashed line) and 
~ —2. 3422 Tjd 2 for = 0.05 (solid line), in excellent agreement with 
numerical values obtained from equation f[22l . 



Using the saddle-point approximation and Stirling's formula (|15p . one can 
prove that 



lim g 

i-^oo \a d 



lr 



X 



e -(lo/d 



see [HI We now apply this result to equation fl25|) and arrive at 

,2 1 



( d 2 \ 
lim exp ( —V— I = (( 

£^oo F \ T J V 



V/T _ (JJ/T _ A \ p -Qr/d 



1)e -Qr/dy_ 



+ c 



(26) 



(27) 



Hence, the asymptotic value of the potential V(£) is given by the following 
exact expression: 



V. = -lln(l 



T 



,I//T _ (e !7/T _ j) 



exp 



d 1 



:28i 



For non-adhesive particles with U = or l r = 0, the asymptotic value of the 
potential V(£) is Voo = —(Tjd 2 ) ln(l — 0). For small bulk volume fractions 
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(f) -C 1 of the particles and large binding energy U with e u ^ T 1, we obtain 

V^-2^1n(l + 0e^ T ^). (29) 

We will use equations (128|) and (129]) in section [5] to calculate the effective 
adhesion energy of surfaces. In the following section HJ we determine the 
global minimum of the potential V(£). 



4 Global minimum of the surface interaction 
potential 

In the present calculation we have chosen the free energy reference state 
in such a way that the effective surface interaction potential vanishes at 
surface contact £ = 0, i.e. V(£ = 0) = 0, see equation ([7]). For surface 
separations < £ < d, the potential V(£) increases linearly with £, i.e. V(£) = 
£T(j)/[(l — 4>)d 3 ], since F gc = for these separations. The potential V(£) 
decreases for separations d < £ < d+l r , and increases again for d+l r < £ < 2d. 
The potential thus attains a minimum at £ — d+l r , see figures [5] and [31 The 
value Vmin = V(£ = l r + d) at this minimum can be calculated again from 
equation (122]) . For £ = d + l r , equation (flOj) reduces to Z\ = e 2U ^ T l r /A and 
Z n = for n > 2 since only a single particle fits into the column. Insertion 
of these results into equation (]2"2l) leads to 



VU. = ~ln 



L (b ( 2U 

1 + "77^ ex P -7^- + 



dl-6 \ T 1 - 



+ 5^ + 1)t^-T- (30) 



d 2 V d 7 1- 



The potential V(£) thus has two local minima, one located at £ = with 
V(£ = 0) = and the other at £ = l r + d with V(£ = l r + d) = V min . The 
two minima result from the interplay of depletion interactions and adhesive 
interactions. The global minimum of V(£) is located at £ = d + l r if < 0, 
i.e. for 

e^ + ^Ce 2l//T -e^/ d >0 (31) 

with ( = 0/(1 — 0), see equation (123]) . The inequality (]3T|) is fulfilled for 
sufficiently large particle binding energies U. For small binding energies U, 
in contrast, the potential V(£) has its global minimum at surface contact 
£ = 0, see figure [31 

In the experimentally relevant case of small particle bulk volume fractions 
<p 1 and large binding energy U with e u ^ T 3> 1, equation (l30j) reduces to 

K^-^lnfl + ^e 2 ^) (32) 
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Figure 3: Rescaled surface interaction potential Vd 2 /T, given by equation 
fl22|) . versus the rescaled surface separation £/d. The bulk volume fraction of 
the particles here is <fi = 0.1, the binding range is l T = 0.3 d, and the binding 
energy is U — 0.9 T (dashed line) and U = 0.6 T (solid line). The global 
minimum of the potential V(£) is located at the separation £ = d + l r for 
U = 0.9 T (dashed line) and at surface contact £ — for U = 0.6 T (solid 
line) . 



and the inequality f[3"Tj) simplifies to 

a>fm( 1 + £). (33) 

and, thus, to a relation that is independent of the particle bulk volume frac- 
tion 0. 

5 Adhesion energy 

In this section, we assume that the binding energy U of the particles is 
sufficiently large so that the inequality (I3TT) is fulfilled. The global minimum 
of the interaction potential V(£) then is located at £ = d + l r . The minimum 
value V min = V{£ = d + l r ) is given by equation (I3"U|) . For large surface 
separations £ ^> d, the potential V(£) attains a constant value given 
by equation fl28]) . The difference between the asymptotic and the minimum 
value of the potential V(£) is the effective adhesion energy 

W = V OQ - V min (34) 




U = 0.6T 



U = 0.9T 
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of the surfaces. The effective adhesion energy is the minimal work that has to 
performed to bring the two surfaces far apart from the separation i — d + l r . 
From equations (1281) . (130]) . and (134|) . we obtain the exact result 



T 

W = ^ ln 

or 



1 + 



d 1 - 



'2U 

,<>X1,) 1 ^r + r 



d 2 w 



1 - 



,!7/T _ (( , K /T _ 



cxp 



d 1 



(35) 



In figure HI the adhesion energy W is plotted as a function of the particle bulk 
volume fraction 0. Interestingly, the adhesion energy W exhibits a maximum 
at an optimal bulk volume fraction 0* of the particles. 
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Figure 4: Rescaled adhesion energy Wd 2 /T as a function of the particle bulk 
volume fraction 0. The solid line corresponds to the exact result (|35|) . and 
the dashed line to the approximation (1361) . The particle binding energy here 
is U = 5T and the binding range is l r = 0.3d. The adhesion energy has a 
maximum at = 0* with 0* »s e~ u ^ T d/l r . 



For small bulk volume fractions 1 and large binding energy U with 
e u/T ^ aS y m ptotic value and minimum value of V(^) are approxi- 

mately given by equations fl29l) and fl32l) . respectively The adhesion energy 
W = Vqo — Vmin then simplifies to 



T l + 0e 2 ^ T / r /rf 
d 2 11 (i + e ^/ r /(i) ; 



rr ^, :, 2 - ( 36 ) 



This expression is identical with our previous result obtained from a virial 
expansion in up to second order terms [8]. For <C 1 and e 17 ^ ^> 1, the 
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adhesion energy (1561) is a good approximation of the exact result (|3"oT) . see 
figures H] and From equation (I3T)|) . we obtain the approximate expression 

f ^^"^ (37) 

for the optimum bulk volume fraction <ft* at which the adhesion energy W 
attains its maximum value. 

The adhesion energy f[3"oT) can be understood as the difference of two Lang- 
muir adsorption free energies per fluid column, or pair of apposing binding 
sites [8j: (i) the adsorption free energy (T/d 2 ) In (l + q <p e 2U ^ T ) for small 
surface separations at which a particle binds both surfaces with total bind- 
ing energy 2U, and (ii) the adsorption free energy (T/d 2 ) In (l + q(pe u ^ T ) 
for large surface separations, counted twice in (I3T)|) because we have two sur- 
faces. These Langmuir adsorption free energies result from a simple two-state 
model in which a particle is either absent (Boltzmann weight 1) or present 
(Boltzmann weights q<pe 2U ^ T and q(pe u ^ T , respectively) at a given binding 
site, see e.g. (2D]- The factor q depends on the degrees of freedom of a single 
adsorbed particle. In our model, we obtain q = l r /d. 

To assess the quality of approximate expression (136|) . we analyze its rela- 
tive error in reference to the exact result (|35l) . The relative error is the mag- 
nitude of the difference between the exact result (1351) and the approximate 
expression (1361) divided by the magnitude of the exact result (1351) . Figure [5] 
shows parameter regions in which the relative error of the expression (I3^|) 
is smaller or larger than 1%, 2% and 5%, respectively. In this example, the 
binding range is l r = 0.3. For intermediate and large binding energies with 
U > 6 T, we find that the relative error of the approximate expression fl36|) 
is smaller than 1% in a broad range of volume fractions </>. 



6 Potential barrier 



For large binding energies U with e u ^ T ^> 1, the effective interaction potential 
has a barrier at surface separations d + 2l r < £ < 2(d + l r ), see figure [2j At 
these separations, only a single particle fits between the surfaces, but this 
particle can just bind one of the surfaces. The particle thus 'blocks' the 
binding site at the apposing surface. 

The potential barrier attains its maximum value Vb a = V(£ = 2d) at the 
separation £ = 2d, see figure El From equation (1231) . we obtain 



ba 



1 + 



(4 ( ^_ 1)+1 )_i_ exp (_i_) 



+ 2 



T 

d 2 ! 



(38) 
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Figure 5: Relative error of the approximate expression (136|) for the binding 
range l r = 0.3/, d of the particles. In the parameter region above the dotted 
line, the relative error is smaller than 1%. Below this line, the relative error 
is larger than 1%. The relative error is smaller than 2% above the dashed 
line, and smaller than 5% above the solid line. For binding energies U > 6T, 
the simple expression (1361) approximates the exact result fl35l) very well since 
the relative error is smaller than 1% for a broad range of volume fractions (p. 



For <p 1 and e u ^ T ^> 1, we get 

^4>„(l + 2^t) (39 ) 

The barrier height = \4>a ~ Voo then is 

T {l + (j ) e u / T l T /d) 2 , 



since the asymptotic value is given by equation (129]) in this limiting case. 
The width of the barrier is approximately Zba ~ d, see figure [2J Equation PU]) 
is again identical with our previous result obtained from a virial expansion 
in (f) up to second order terms [S]. 
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7 Binding probability 



Another quantity of interest here is the binding probability 




defined as the probability that the separation of the closest particle from 
a column base is smaller than the binding range l r . In other words, the 
binding probability n s is the probability of finding a particle bound to one of 
the bases. The binding probability corresponds to the surface coverage in the 
case of a three-dimensional gas of particles between two parallel attractive 
surfaces. 

Equations (JTJ) - (J5j) imply that the binding probability can be calculated 
by differentiation of the grand potential F gc with respect to binding energy U, 

i.e. n s = — | Since the grand potential density f gc given by equation 

(J2TJ) does not depend on the binding energy U, the binding probability can 
also be obtained from the effective surface interaction potential via 



With the exact expression ([25]) for the interaction potential V, the binding 
probability n s can be determined numerically for any finite separation £. 
In figure El the binding probability n s is plotted as a function of surface 
separation £ for three different volume fractions around the optimal volume 
fraction 0* at which the adhesion energy W is maximal. In the vicinity of 
0*, the binding probability at large separations £ > 2(d + l r ) is sensitive to 
small variations of 0, while the binding probability at separations £ in the 
surface binding range d < £ < d + 2l r remains practically constant at almost 



For small bulk volume fractions of the particles and large particle bind- 
ing energies U, the asymptotic and minimum value of the interaction po- 
tential V(£) are given by equations fl29|) and fl32|) . respectively. From these 
equations and relation (HS]) . we obtain the approximate expressions 



for the binding probability at large surface separation £ and 
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Figure 6: Binding probability n s , calculated numerically from equations (j42p 
and fl25l) . as a function of rescaled surface separation Ijd where d is the 
diameter of the adhesive particles. The binding energy here is U — 7T, 
the binding range is l r = 0.3 d and the particle bulk volume fraction is = 
0.001 < 0* (dashed line), = 0.003 » 0* (solid line) and = 0.007 > 0* 
(dotted line). The optimal volume fraction 0* at which the adhesion energy 
becomes maximal is given by equation (1371) . 



for the particle binding probability at the binding separation £ = d + l r of the 
surfaces, with the optimum bulk volume fraction 0* given in equation fl37j) . 
These expressions correspond to the well-known Langmuir adsorption equa- 
tion [20J. At the optimal volume fraction 0*, the particle binding probability 
for unbound and bound surfaces is n Sj00 = 1/2 and n Sjmin rs 1, respectively. 
Bringing the surfaces from large separations £ > 2(d + l r ) within binding 
separations d < £ < d + 21 r thus does not require desorption or adsorption of 
particles at = 0*. 



8 Conclusions 

We have considered one-dimensional gas of hard-sphere particles with at- 
tractive boundaries, a novel extension of the Tonks model [Ij. We have 
solved this model analytically in the whole range of parameters by explicit 
integration over the particles' degrees of freedom in the partition function. 
In contrast to other studies on one-dimensional models for hard spheres 
HH HH [201 EH [221 I231 [2H [251 126], we have focused on the boundary con- 
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tribution to the free energy of the system, which corresponds to the effective, 
particle-mediated interaction potential between the boundaries, or surfaces, 
see figures [2] and [3j The effective adhesion energy obtained from the inter- 
action potential depends non-monotonically on the volume fraction of the 
particles in the bulk, see figure HI The adhesion energy exhibits a maximum 
at an optimum volume fraction, which can lead to reentrant transitions in 
which the surfaces first bind with increasing volume fraction 0, and unbind 
again when the volume fraction <fi is increased beyond its optimum value. 

A lattice of such one-dimensional gas columns represents a discrete ap- 
proximation of a three-dimensional gas of particles between two adsorbing 
surfaces, see figure [1] and reference [8]. For small volume fractions <fi and 



short-ranged particle-surface interactions considered here, the gas of parti- 
cles between two well-separated surfaces is as dilute as in the bulk, except for 
the single adsorption layers of particles at the surfaces. At larger volume frac- 
tions, three-dimensional packing effects become relevant. These effects are 
not captured correctly in the one- dimensional model. However, it has been 
pointed out [22] that approximations based on one-dimensional models do 
well in comparison to density functional theories for three-dimensional hard 
sphere fluids confined in planar, non-adsorbing pores [27] . In principle, the 
quality of the one-dimensional approximation can be tested by Monte Carlo 
or Molecular Dynamics simulations, which have been used to study various 
three-dimensional systems of hard spheres confined between non-adsorbing 
surfaces [2B3 M, M, EE E2] • 

For simplicity, and for consistency with our previous publication [8J, we 
have considered here a square lattice of columns between the surfaces. In 
particular, the factor d 2 in the denominator of equation ([7]) is the column 
base area in the square lattice. For a hexagonal lattice of columns, the 
corresponding area is (-\/3/2) d 2 , and the corresponding effective interaction 
potential of the surfaces is thus obtained by multiplying the right hand side 
of equation fT2"2l with a factor 2/v^3 ~ 1.1547. The adhesion energy of 
the surfaces has to be rescaled with the same factor in the case of hexagonal 
lattice of columns, but its functional dependence on the bulk volume fraction 
<f), binding energy U, and binding range l r remains unchanged. 

We have considered an equilibrium situation in which the particles ex- 
change with a bulk solution. For polymers between surfaces, such an equi- 
librium has been termed 'full equilibrium'. In a 'restricted equilibrium', in 
contrast, the polymers are trapped between the surfaces (33J, EU [35], which 
is less likely for the spherical particles considered here. 
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A Canonical partition function 

In this section we calculate the n-particle partition function Z n as given by 
equation @ for I > nd. First, if one notices that 

e ue(i r - yi )/T = 1 + ( e u/T _ ^ Q(/r _ yi) (45) 

and 

e Ue(y n -i+nd + l r )/T = x + ( e U/T _ ^ _ g + ^ + ^ (4g) 

the integral flU]) can be written as a sum of four terms 
with 

ri-nd py n ny 3 ny 2 

h = dy n dy n ^ ... dy 2 dy u (48) 



JO JO JO JO 

t-nd i>y n i>y 3 i>y 2 

dy n Q (y n - i + nd + l r ) / dy n -i . . . / dy 2 / dy x , (49) 
Jo Jo Jo 

rl-nd ry n ny 3 py 2 

h= dy n / dy n _i . . . / dy 2 / dj/i6 (Z r - y x ) , (50) 

JO JO Jo Jo 

'-nd py„ ry 2 

dy n Q (y n - £ + nd + l r ) dyn-!... dy x Q (l r - y x ) . (51) 



The first integral 



I rl—nd | 



and second integral 

j fl—nd 



(n-1)! 



= —((€- nd) n -(£-nd- IS 0(£-nd- U) (53) 
n! 

can be easily calculated. 

To calculate the third integral, we start from 

/ 6 (l r - y x ) dy x = min [y 2 , l r \ = yi~ (?/2 - U) © (?/2 - • (54) 

</ o 
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In the next steps, we find 

2/3 11 

{Vl - (V2 ~ lr) @ (2/2 - l r )) d V2 = 7^3 ~ j ^ 3 ~ ^ 9 (#3 ~~ ( 55 ) 



and 



//1 



- ~ ( '/:; - /<•)"' B ( //.-. - /,. ) ) dy 3 = - - (y 4 - Z r ) 3 9 (y 4 - / r ) 



Iterating these results leads to 



(56) 



1 



i-nd 



° (n- 1)1 J 
The integral J3 can now be evaluated as 



K" 1 - fan - lr) 6 (j/ n - Z r )] dy n . (57) 



h = —At- nd) n - — / (y n - Z,.)™" 1 6 (y n - l r ) dy n 

= —At- nd) n - — {i-nd- l r ) n eU-nd- l r ) . (5£ 

Note that I 2 = h- 

The fourth integral, I4, can be brought to the form 



j pi—nd 



(n-l)\Jo 



Q{y n -£ + nd + l r ) [y*- 1 - (y n - lr) 11 ' 1 6 (y n - l r )] dy r , 



(59) 

if one uses again (|5I]) and iterates the integration as in (1531) and (f5B"j) . The first 
term on the right hand side of equation (1591 is equal to J 2 , see equation (1551) . 
Thus I4 — I2 — h, where 

© (l/n - £ + nd + l r ) 6 (y n - l r ) (y n - lr) 11 ' 1 dy n (60) 



° (n-l)!j 

To determine the integral -Z5, one has to distinguish three cases: (i) for £ > 
nd + 2Z r , we obtain 

7 5 = 1 {{£ -nd- l r ) n -{i-nd- 2l r ) n ) , (61) 
n! 

(ii) for nd + l r < £ < nd + 2Z r , we obtain 

h = —Al-nd- l r ) n (62) 
n! 
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and (iii) for I < nd + l r one gets I5 = 0. In summary 

1 



rv. 



— nd — l r ) n (£ — nd — l r 



+— (£-nd- 2l r ) n G(£-nd- 2l r 



nl 



(63) 



If one now gathers the results ([52]) . ([55]) . (j581) . ([55]) and returns to equation 
(I47p . one obtains the partition function 

Z « = 7^7 [ J i + 2 ^ U ' T ~ X )( J i " J 2) + (^ /T - l)Vi - 2J 2 + J 8 )] (64) 



with 



J x = (i- nd) n 9 {£ - nd) , 
J 2 = (£ - nd - l r ) n e(£-nd- l r ) 
j 3 = (£-nd- 2l r ) n Q(£-nd- 2l r 
This result can be written as 
1 



A n n\ 



(e u/T - 1) {£-nd- 2l r ) n Q(£-nd- 2l r ) 
2e u/T (e u/T -l){£-nd- l r ) n Q (£ - nd - l r ) 
+e 2U/T {£ - nd) n Q{£- nd) 



which simplifies to equation (jTTJ 



(65) 
(66) 
(67) 



(68) 



B Proof of equality (1261) 



Here, we explore the asymptotics of the function g(£/d,lo/d,() defined in 
equation f[24l) and, hence, prove the equality ([261) . To simplify the notation, 
let £/d = N, where N is a large integer number, and l$/d = A. Then 



N 



An 

g (jV, A, C) = V ±-e~ {N - nK (N - n - X) n 9 (N - n - A) . (69) 
f n! 



n=l 



In the next step, we introduce the auxiliary function 



x{ n ) — —{N — n)( + n ln( + n In (AT — n — A) — nlnn + n ln(27m) (70) 
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and rewrite the function g(N, A, () given by equation ([611 in the form 



N 



g = ^e^e(iV-n-A), (71) 



n=l 



using Stirling's formula (j!5p . For large TV, one can replace the sum on the 
right hand side of equation (I7TT) by an integral and write 



AT— A 

e x(n) dn. (72) 

The function x( n ) has a global maximum at n = uq with 

"o « (iV - A) ^ (73) 

for large JV. Note that for large numbers iV, the location n « (iV — A) of 
the global minimum scales linearly with N. If we now expand the function 
x(n) around the point n up to second order terms and apply the saddle-point 
approximation, we get 



g m " / exp 



-x"(n ) (n-n ) 2 



dn. (74) 



A simple change of variables m = n — n$ leads to 



y(JV-A)(l-*) „ 

g « e x(no) / e x (no)m /2 dm (75) 

J-{N-X)cf>+1 

where we have used no ~ (N — A) <f), which follows from equation (175]) and 
relation (l25j) between variables ( and 0. In the limit of large N, we thus 
obtain 

/oo „ 
e X (n )m 2 /2 dm (yg) 
-oo 

because < (ft < 1. Now, we can calculate the Gaussian integral in (176]) to 

get 

From the definition (170]) of the auxiliary function x(n) and equation (175]) for 
the point n = no at which function xi n ) nas its global maximum, we get 

xK)~-AC-^hi(27m ) (78) 
21 



and 

X(n ) - (79) 
no 

Note that x" ( n o) < 0, and that the function x( n ) nas indeed a maximum at 
n = n . Combining equations (J77|) . (1751) and (1791 leads to 



e -AC 



° l + C 

for large N values and, thus, to equation (|26i) quod erat demonstrandum. 



,0) 
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